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By evaluating multipole susceptibility of a seven-orbital impurity Anderson model with the use of 
a numerical renormalization group method, we discuss possible multipole states of actinide dioxides 
at low temperatures. In particular, here we point out a possible scenario for multipole ordering in 
americium dioxide. For Am 4+ ion with five 5/ electrons, it is considered that the ground state is F^ 
doublet and the first excited state is F% quartet, but we remark that the / ground state is easily 
converted due to the competition between spin-orbit coupling and Coulomb interactions. Then, we 
find that the F% quartet can be the ground state of AmC>2 even for the same crystalline electric 
field potential. In the case of quartet ground state, the numerical results suggest that high-order 
multipoles such as quadrupole and octupole can be relevant to Am02. 



PACS numbers: 75.20.Hr, 75.40.Cx, 71.70.Ch, 71.27. +a 



I. INTRODUCTION 

Actinide dioxides with the fluorite structure of the 
space group Fm3m have been studied intensively for 
more than fifty years both from experimental and the- 
oretical sides.— «2 A typical target material is UO2 mainly 
due to its technological importance as a nuclear reactor 
fuel and a heterogeneous catalyst. As for theoretical re- 
search, a clear picture of the electronic structure in UO2 
has been obtained. In fact, neutron scattering results 
show that UO2 is a noncollinear antiferromagnet below 
30. 8K.- Detailed analysis of core photoemission spectra 
has suggested that UO2 is an insulator of Mott-Hubbard 
type4 The crystalline electric field (CEF) states of UO2 
have been also determined. 5 

Concerning NpC-2, over fifty years since 1953^ it has 
been known to exhibit a mysterious low-temperature or- 
dered phased Several phenomenological works on the 
ordered phase have claimed a key role of octupole de- 
gree of freedom.2ii2.iI.I2.I2 The CEF states of Np0 2 have 
been determined by neutron scattering experiment^ Re- 
cently, the octupole order has been strongly supported by 
17 0-NMR experimentA 5 - and by inelastic neutron scatter- 
ing study.i2 In order to understand why such high-order 
multipole ordering appears, it is necessary to proceed 
to the research in a microscopic level. For this issue, 
it has been shown that octupole order actually charac- 
terizes the ground state of NpC>2 by the analysis of an 
/-electron model on the basis of a j-j coupling scheme 
on an fee lattice i 17 ' 18 

On the other hand, PuC>2 is known to be a semicon- 
ductor with magnetic susceptibility which is almost in- 
dependent of temperature up to 1000 K, since the CEF 
ground state is Tf singlet and the first excited state is 
triplet with the large excitation energy as 123 meV i 19 ' 20 
Thus, from the viewpoint of magnetism, PuC>2 did not 
attract much attention. 

Now let us turn our attention to Am02- Probably due 
to the difficulty in the treatment of this material with 
high radioactivity, we cannot find lots of experimental 



results on Am02. In 1969, Mossbauer isomer shift in 
Am02 was measured^ After that, the magnetic suscep- 
tibility was measured and the peak was found around at 
15 K. 22 Naively thinking, such a peak seems to suggest 
the signal of antifcrromagnetic ordering, while neutron 
diffraction measurement did not detect antifcrromagnetic 
order in agreement with the Mossbauer measurement^ 
This situation looks similar to that of NpC>2. Namely, 
multipole degree of freedom seems to be a key issue in 
A111O2 to reconcile experimental results, as has been pro- 
posed for NpC-2 in the context of octupole ordering. 

Here we note that the CEF ground state of A111O2 
was considered to be Ty doublet from the experimen- 
tal resultS ) 22 ' 24 i 25 but it could not bring higher multi- 
poles. This is in sharp contrast to the case of NpC-2 with 
the confirmed CEF ground state of T% quartet. Due to 
the CEF analysis of actinide dioxides, the doublet 
ground state has been suggested^ Thus, it seems to be 
the mainstream in the research of actinide dioxides to 
clarify a mechanism which explains the disappearance of 
antifcrromagnetic order for the Ty doublet ground state. 

However, we believe that there still exists an al- 
ternative scenario on the basis of multipole ordering 
in AmC>2, when we recall the fact that the CEF 
ground state of the / 5 electron system is easily con- 
verted due to the competition between spin-orbit cou- 
pling and Coulomb interaction!^ 7 . This point has been 
also discussed by the present author to propose a 
possible scenario which explains the change of the 
CEF ground state among Sm-based filled skutterudite 
compounds, 28 since trivalent Sm ion includes five 4/ 
electrons. It has been experimentall y 29 ' 30 ' 31 ' 32 and 
theoreticall y 29 ' 30 ' 31 ' 32 ' 33 shown that the angular momen- 
tum coupling of the 5/ states of Am is situated between 
the LS and j-j coupling limits for many chemical sit- 
uations, albeit closer to the j-j limit. Thus, a strong 
spin-orbit coupling is present in the 5/ states of / 5 and 
/ 6 configurations of Am that can compete with Coulomb 
interactions. Accordingly, it is believed to be meaning- 
ful to pursue a possibility of multipole ordering in A1T1O2 
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with the T% quartet ground state. 

In this paper, we show that multipole order is possi- 
ble in Am02, when we appropriately take into account 
both spin-orbit coupling and Coulomb interaction in the 
/-electron terms for the CEF ground state. In the case of 
five / electrons such as Am 4+ ion, the ground state is eas- 
ily converted between and T%, when spin-orbit cou- 
pling and Coulomb interaction compete with each other. 
Then, the quartet can be the ground state, even if 
the CEF potential is unchanged. In order to see what 
type of multipole is relevant, we evaluate the multipole 
susceptibility of the Anderson model by using a numeri- 
cal renormalization group method. We find that higher- 
order multipoles are actually relevant to Am02 within 
the present calculation. 

The organization of this paper is as follows. In Sec. II, 
we discuss the local /-electron state emerging from the 
competition among the Coulomb interaction and spin- 
orbit coupling under the CEF potential. In particular, 
the change of the / 5 ground state is explained in detail. 
In Sec. Ill, we show the model Hamiltonian. In order to 
discuss multipole properties, it is necessary to define the 
multipole operator. Here we explain the description of 
the multipole as spin-charge complex one electron oper- 
ator. Then, we briefly explain the numerical technique 
used in this paper. In Sec. IV, we show the results of 
the multipole state for the case of n = 2 ~ 5, where n is 
the local /-electron number. In particular, the results for 
n=5 are discussed in detail. Finally, the paper is sum- 
marized in Sec. V. Throughout this paper, we use such 
units as fh=kB=l. 

II. LOCAL /-ELECTRON STATE 

Let us first discuss the local /-electron states of ac- 
tinide ions. The local Hamiltonian should be composed 
of three parts as 

H\ oc = H so + Hi n t + HcEF- (1) 

The first term denotes the spin-orbit coupling, given by 

H S o A ^ ' ^ ' Cm,a-,m', a' fmtr fra'a' ; (2) 
m,m' tr,cr' 

where <r=+l (—1) for up (down) spin, f ma is the an- 
nihilation operator for / electron with spin a and z- 
component m of angular momentum £=3, A is the spin- 
or bit coupling, C m,±i,m,±i = ±m/2, Cm±i,=Fi,m,±i = 
y/12 — m(m ± l)/2, and zero for the other cases. 

The second term indicates the Coulomb interaction 
among / electrons, expressed as 

Hint = ^ ' ^ 1 1mim2,m3m,i fmicrfmya' fmstr' fmjai (3) 
mi ~m4 a,a' 

where the Coulomb integral Irn 1 ,m 2 ,m 3 ,m 4 is given by 

6 

Imi,m 2 ,m 3 ,m4 ^^F k c k {mi,m i )c k {m 2 ,m 3 ). (4) 

fc=0 



Here F k is the radial integral for the fc-th partial wave, 
called Slater integral or Slater-Condon parameter—^ 
and c k is the Gaunt coefficient i 36 ! 37 Note that the sum 
is limited by the Wigner-Eckart theorem to even values 
(fc=0, 2, 4, and 6). 

The third term is the CEF potential, given in the one- 
electron potential form as 

H C EF = ^5 m ,m'/L/m'f, (5) 

m,m' cr 

where B m>m i is the CEF potential. Since the fluorite 
structure belongs to Oh point group, B m , m i is given by 
using a couple of CEF parameters B\ and Bq for angular 
momentum 1=3 as 3 ^ 

B 3 . 3 = B_ 3 ._3 = 180B 4 ° + 1805°, 
B2.2 = B-2,-2 = -420B2 - 1080Bg, 

= = 60B 4 ° + 27005°, 

B .o = 360B1 - 3600Bg, W 
Sa!-i - £-3.1 = 6QVl5(B2 - 215°), 
£2,-2 = 300^2 + 7560Bg. 

Note the relation of B m ^ m i—B m i m . Following the tradi- 
tional notation^ we define 

B° 4 = Wx/F(i), B° 6 = W(l - M)/F(6), (7) 

where x and the sign of W specify the CEF energy 
scheme, while \W\ determines the energy scale for the 
CEF potential. Concerning non-dimensional parameters, 
F(4) and F{6), we use F(4)=15 and F(6)=180 for £=3. 38 

Here we briefly explain the parameters of the local 
Hamiltonian. Concerning Slater-Condon parameters, 
first we set F°=10 eV by hand, since we are not interested 
in the determination of the absolute value of the ground 
state energy. Others arc determined so as to reproduce 
excitation spectra of U 4+ ion with two 5/ electrons^ 
Here we show only the results: F 2 =6.36 eV, _F 4 =5.63 
eV, and F 6 =4.13 eV. As for spin-orbit coupling A, we 
use the values of actinide atoms such as A=0.235 eV (U), 
0.272 eV (Np), 0.311 eV (Pu), and 0.351 eV (Am). 

For the estimation of the CEF parameters, let us sum- 
marize the CEF energy levels of actinide dioxides. For 
UO2, the ground state is triplet and the first excited 
state is Tg doublet with the excitation energy 150 meV£ 
For Np02, the ground and first excited states are, re- 
spectively, r g ^ and r g ^ quartets with the excitation 
energy 55 meV For PuC-2, the ground state is sin- 
glet, while the first excited state is T 4 triplet with the 
excitation energy 123 meVi 19 ' 20 

Now we set the CEF parameters for actinide dioxides in 
the present notation. For the purpose, first we estimate 
W and x so as to reproduce the CEF scheme of UO2 
by the / 2 electrons state. This is not a difficult task, 
since the CEF parameters are easily restricted from the 
experimental results. After that, among the values of 
W and x appropriate for UO2, we further restrict the 
values of W and x which can also reproduce the results 
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FIG. 1: CEF energy level schemes for tetravalent actinide 
ions. Parameters used here are explained in the maintext. 



for Np0 2 and Pu0 2 . Note that the CEF states for f 3 and 
/ states are almost reproduced by using the parameters 
of f 2 electron state. Since the CEF term is just given 
by the one-electron potential, the CEF effect is not so 
drastically changed among the materials with the same 
crystal structure, even though the /-electron number is 
different. 

After some calculations, we determine W=— 10.5 meV 
and £=0.62. The results for the CEF level schemes are 
summarized in Fig. 1. It is observed that the CEF states 
of UO2, Np02, and PuC-2 are well reproduced. Note here 
that for the same CEF parameters, the ground state for 
is found to be doublet and the excited state is 



Am 44 

quartet with the excitation energy of about 50 meV. 
This is consistent with the previous theoretical results 
on the CEF states of Am02 obtained by more detailed 
calculations^ 

Naively thinking, there occurs ordering of magnetic 
moment originating from ground state for A111O2, but 
it seems to be contradict with neutron diffraction study. 
In order to resolve such contradiction, there are two ways: 
One is to consider a mechanism which explains the dis- 
appearance of magnetic moment even in the Ty ground 
state. Another way is to reconsider the local /-electron 
term by focusing on spin-orbit coupling and Coulomb in- 
teraction. Here we propose an alternative scenario on the 
basis of the second direction. 

Thus far, we have simply assumed that Coulomb inter- 
action is not changed among different actinide ions, but 
in actuality, they may be changed. Concerning spin-orbit 
coupling, we have used the values in actinide atom, but it 
may be also changed in the ionic states. Then, we point 
out that for / 5 systems, the CEF states are sensitively 
changed by the competition between Coulomb interac- 
tion and spin-orbit coupling in comparison with other 
values of local /-electron numberj 27 i 28 

In order to understand that the effect of Coulomb in- 
teraction and spin-orbit coupling appears in the CEF pa- 
rameters, here we express the CEF parameter B\ by us- 
ing the so-called Stevens factor as 




FIG. 2: The CEF energy for n=5 vs. r. The meaning of r is 
explained in the maintext. 



where Ak is the parameter depending on materials, (r k ) 
denotes the radial average of local /-electron wavefunc- 
tion, n denotes the local /-electron number, J is the 
total angular momentum of the ground state multiplet, 
and /3j indicates the Stevens factor, which is one of co- 
efficients appearing in the method of Stevens' operator 
equivalent^ 

For the case of n=5, it is well known that the 
ground state multiplet is characterized by J— 5/2. Af- 
ter lengthy calculations, for n=5 and J=5/2, we can ob- 
tain that f3^ 2 =(13/ 21) f3 3 in the LS coupling scheme and 



/3g/2=™(ll/7)/33 in the j-j coupling scheme^! where (3$ 
denotes the Stevens factor for £=3, given by /3 3 =2/495. 
It should be noted that the sign of f3j^ is changed be- 
tween the LS and j-j coupling schemes, suggesting that 
the ground state is converted, when Coulomb interaction 
and/or spin-orbit coupling are changed. 

The modification of Coulomb interaction and spin- 
orbit coupling is closely related to the picture for multi- 
/-electron state. The CEF level schemes in Fig. 1 is qual- 
itatively understood by a j-j coupling scheme. Namely, 
by assuming that the effective Hund's rule coupling is 
smaller than the CEF level splitting, we simply accom- 
modate plural numbers of / electrons in the levels of 
r^" ground and Tj excited states. Then, we can easily 
reproduce all the CEF level schemes of tetravalent ac- 
tinide ions. However, the conversion of the CEF ground 
state for n=5 indicates that the actual situation should 
be slightly shifted to the side of the LS coupling scheme. 
We note that due to such a shift, the CEF states for n=2, 
3, and 4 are not qualitatively changed, while in the case 
of rt=5, the ground state is converted. 

Here we emphasize that such a ground-state conversion 
occurs in the region of realistic values of Coulomb inter- 
action and spin-orbit coupling, as has been pointed out 
in the discussion of the CEF states of Sm-based filled 
skutterudites^ In fact, as mentioned in Sec. I, it has 
been shown that the angular momentum coupling of the 
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5/ states of Am is situated between the LS and j-j 
coupling limits, but it is rather closer to the j-j cou- 
pling limit . 29 i 30 ' 31 i 32 ' 33 Thus, it seems to be reasonable 
to change slightly spin-orbit coupling and/or Coulomb 
interaction in the case of tetravalent Am ion. 

If we resort to first-principles calculation, we may de- 
termine correctly F-? and A, but it is out of the scope of 
the present paper. Here we simply introduce an artificial 
parameter r to control Coulomb interactions and spin- 
orbit coupling as F k — > rF k and A — > A/r with r > 1. 

In Fig. 2, we depict the energies of H\ oc as functions of r 
for n=5. Wc find that the conversion of the ground state 
occurs around at r w 1.3. Note that in other values of n, 
the ground states are not changed. Thus, we can obtain 
the r^" quartet ground state for Am 4+ . In the following 
sections, we will discuss how the multipole states actually 
appear when we change the value of r for Am ion. 



III. MODEL AND METHOD 

A. Anderson model 

Now we include the hybridization between localized 
and conduction electrons. The Hamiltonian is the An- 
derson model, given by 

H =E £ * c L c ^ + E ( v ™ c Lf™° + kc -) + ff io C) (9) 

k.cr k,cr,m 

where £k denotes conduction electron dispersion, Cko- in- 
dicates the annihilation operator for conduction electron 
with momentum k and spin a, and V m is the hybridiza- 
tion between conduction and / electrons. 

Note that we consider only a u single conduction band 
with xyz symmetry composed of oxygen 2p electrons. 
Since oxygen ions surrounding actinide ions are located 
in the [1, 1, 1] direction, there should exist a conduc- 
tion band composed of 2p electrons with xyz symmetry. 
This picture seems to be consistent with band-structure 
calculation^ but the ignorance of t\ u and t2 U bands is 
just assumption. Here wc note that the hybridization 
occurs between the states with the same symmetry of lo- 
cal /-electron state. Since the a u conduction band has 
xyz symmetry, we set V2=—V-2=V and zero for other m. 
Hereafter, a half of the bandwidth of a u conduction band 
D is set as the energy unit, i.e., D=l eV. We fix V as 
V/ D=0.05 throughout this paper. Note that in order to 
adjust the local /-electron number n, we appropriately 
change the chemical potential in the actual calculation, 
although we do not explicitly show such a term. 



B. Multipole operator 

In order to discuss multipole properties, it is nec- 
essary to define the multipole operator X for / 



electroni 43 i 44 i 45 i 46 i 47 i 48 In general, X is expressed as 

X = 5>( fc )fW (io) 

where k is a rank of multipole, 7 is a label to express 

* (k) 

Oh irreducible representation, and T~j is cubic tensor 
operator, given by Ty = J2 q G^ k } q Tq k ^ . Here an integer q 
runs between — k and k, Tq k ^ is spherical tensor operator, 

(k) 

and G\.q is the transformation matrix between spherical 
and cubic harmonics. We determine p^ later. 

In order to obtain explicit expression of the sphcri- 

" (k) 

cal tensor operator Tq , it is convenient to convert the 
/-electron basis from (m, a) to where j is the to- 

tal angular momentum and \i is the z-component of j. 
When we define fj^ as the annihilation operator for / 

* (k) 

electron labeled by j and /1, we obtain Tq in the second- 
quantized form as 

Note that there are no components between different val- 
ues of j, since the matrix for total angular momentum 
is block-diagonalized in the basis of (J, fj,). The matrix 

element of T^^, is calculated by the Wigner-Eckart the- 
orem as 

where {j^\jn'kq) indicates the Clebsch-Gordan coeffi- 
cient and (j||T( fc )||j) denotes the reduced matrix ele- 
ment for spherical tensor operator, given by (j\ \ \j)= 
v / (2j + k + l)!/(2j - k)\/2 k . Note that k<2j and the 
highest rank of /-electron multipole is 7. 

(k) 

Let us now determine the coefficient p\ . In order to 
discuss the multipole state, it is necessary to evaluate 
the multipole susceptibility in the linear response theory. 
However, multipoles belonging to the same symmetry are 
mixed in general, even if the rank is different. In addi- 
tion, multipoles are also mixed due to the CEF effect. 

(k) 

Thus, we determine p\ by the normalized eigenstate of 
susceptibility matrix 

= E ._ E . (i\m k) -p^}\3) 

id 3 1 

x (j\[ff-p^]\i), (13) 

where Ei is the eigenenergy for the i-th eigenstate \i) of 
H, T is a temperature, Py ] =Y,i e~ E ^ T {i\f\ k) \i)/Z, and 
Z is the partition function given 

b y z=Y.i e ~ Ei/T - Note 

that the multipole susceptibility is given by the eigen- 
value of the susceptibility matrix. 
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C. Method 

In order to evaluate multipole susceptibility of the An- 
derson model, we employ a numerical renormalization 
group (NRG) method^ in which momentum space is 
logarithmically discretized to include efficiently the con- 
duction electrons near the Fermi energy and the conduc- 
tion electron states are characterized by "shell" labeled 
by N. The shell of N=Q denotes an impurity site de- 
scribed by the local Hamiltonian. 

In the NRG method, we transform the Hamiltonian 
into the recursion form as 

c Na C N+lcr + c\j +la CNa), (14) 



where A denotes a parameter for logarithmic discretiza- 
tion, cjvc indicates the annihilation operator of conduc- 
tion electron in the iV-shell, and tjy is the hopping of 
electron between N- and (N + l)-shells, expressed by 



(1 + A-!)(1-A 



-N-l\ 



(15) 



2y/(l- A- 2Ar -!)(l- A- 2Ar - 3 ) 
The initial term Hq is given by 

H a = A-^IH^ + V(clf ca + f^coa)]. (16) 

a 

Each component of multipole susceptibility matrix 
cq. (fT5|) is evaluated by using the renormalized state. 
Then, the multipole state is defined by the eigen states 
of eq. (|13|) . We note that the temperature T is defined 
as r^A"^- 1 )/ 2 in the NRG calculation, where N is the 
number of the renormalization step. Due to the limita- 
tion of computer resources, we keep only M low-energy 
states. In this paper, we set A=5 and M=3000. 

IV. RESULTS 

Now we discuss the multipole state of the Anderson 
model eq. ©. In Figs. 3(a)-(c), we show Tx 7 vs. T for 
n—2, 3, and 4, where x-y is the eigenvalue of the multipole 
susceptibility. We use the same values of the parameters 
in the Hamiltonian for n=2, 3, and 4. The values of 
the spin-orbit coupling are A=0.235 eV, 0.272 eV, and 
0.311 eV for n=2, 3, and 4, respectively. The eigenstates 
are classified by irreducible representation of Oh point 
group. Here we use short-hand notations such as "3g" 
and "5u" , which denote and , respectively, in the 
Bethe notation. Note that "lu" does not appear among 
multipoles up to rank 7. 

For n—2 (UO2), we find the optimized state is la- 
beled by 4u, which is mainly composed of dipole (about 
92%). The secondary components are quadrupoles (5g 
and 3g), but X5g an d X3g is smaller in one order in com- 
parison with Xiu- For n=3 (NpO^), as easily deduced 
from the Tg quartet ground state, we find varieties of 
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FIG. 3: T Xl vs. temperature T for (a) UO2, (b) Np0 2 , and 
(c) PuQ 2 . 



multipoles. Among them, the primary component is 4u, 
which is mainly composed of dipole (about 96%). The 
secondary one is 2u octupole. In the third group, another 
4u and 5u are almost degenerate. We note that this 4u 
is composed of higher multipole component and the 5u is 
mainly composed of octupole. The fourth component is 
3g quadrupole. We emphasize the existence of octupoles 
(2u and 5u) with significant eigenvalues. 
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is doublet, the component which remain in the low- 
temperature region is 4u, which is composed of dipole 
(about 25%) and octupole (about 75%). We note that 
the octupole component is significantly large in compar- 
ison with the 4u states of UO2 and NpC>2. In any case, 
as we have expected, we find only magnetic 4u moment 
originating from ground state for Am02- 

Then, we increase the value of r in order to move to 
the side of the LS coupling scheme. In Fig. 4(b), we plot 
T\- f of the multipole susceptibility for n=5 and r=1.5 
with the Tg ground state (see Fig. 2). In this case, we 
find that the primary component is 3g, which is mainly 
composed of quadrupole. The secondary component is 
4u and we also find 5g, 4u, and 5u components with 
smaller eigenvalues. We note that some multipoles of 
A111O2 are the same as those in NpC>2, except for 2u and 
5g, although the corresponding eigenvalues are different. 

Now we provide a comment on the value of r, which is 
introduced so as to increase the effect of Coulomb interac- 
tions and decrease the magnitude of spin-orbit coupling. 
Note, however, that the angular momentum coupling of 
the 5/ states of Am is nowhere near as close to the LS 
coupling limit as for Cm . 31 ' 33 Thus, we should not en- 
tirely or strongly suppress the spin-orbit coupling in the 
5/ states of Am. It is acceptable that we change slightly 
spin-orbit coupling and Coulomb interaction. Accord- 
ingly, when the value of r is increased so as to move 
towards the side of the LS coupling scheme, we should 
pay due attention not to suppress spin-orbit coupling too 
much. 



FIG. 4: (a) T\-, vs. temperature T for n=5. The parame- 
ters are the same as those in Fig. 2, except for the spin-orbit 
coupling, (b) T\~t vs. temperature T for n=5 and r=1.5 
with enhanced Coulomb interaction and reduced spin-orbit 
coupling. 



In the present calculation, we cannot determine the 
kind of multipole ordering in actual systems. However, 
the multipoles which remain at low temperatures are the 
candidates which will order in the actual system. In the 
case of Np02, it has been gradually revealed that triple-q 
order of 5u octupole can naturally reconcile several kinds 
of experiments. The 5u octupole is actually included in 
the multipoles in the present calculations, although it is 
not dominant component. 

For n=4 (PUO2), we do not find any significant multi- 
pole component, as easily understood from the singlet 
ground state which is well separated from the magnetic 
excited state. In this sense, from the viewpoint of mag- 
netism, this case does not attract much attention. 

Next we move on to the case of n=5, corresponding 
to Am02- In Fig. 4(a), we show the results for n=5 by 
using the parameters for r=l. Namely, the Coulomb in- 
teractions are the same as those in Figs. 3(a)-(c). The 
spin-orbit coupling is set as A=0.351 eV, which is the 
value for Am atom. In this case, since the ground state 



V. DISCUSSION AND SUMMARY 

In this paper, we have discussed the multipole state 
in the low-temperature region, by analyzing the seven- 
orbital impurity Anderson model with the use of the 
NRG method. We have found the multipole state for 
n=2, 3, and 4, which are not in contradiction to the 
phases observed in UO2, Np02, and Pu02, respectively. 
Note here that we determine the candidates which will 
order in the actual periodic system at low temperatures. 

For the case of n=5, when we use the same parame- 
ters as those for n = 2 ~ 4, we have suggested the phase 
dominated by magnetic moment. However, if we change 
slightly Coulomb interactions and spin-orbit coupling, we 
have found the Tg" quartet ground state for A111O2. In 
this situation, we have shown that the low-temperature 
phase can contain multipoles such as quadrupole and oc- 
tupole. The parameters are artificially introduced here, 
but our purpose is to point out a possibility of the Tg - 
quartet ground state due to the competition between 
Coulomb interactions and spin-orbit coupling. 

Unfortunately, we cannot determine the kind of mul- 
tipole order only from the present calculation, but on 
the basis of the same crystal structure, it is plausible 
that 5u octupole order also appears in A111O2. On the 
other hand, it may be possible to exploit other scenar- 
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ios, e.g., quadrupole ordering, which were invented for 
understanding of Np02- In any cases, the combination 
of phenomcnological theory and microscopic experiment 
will be useful to finalize the kind of multipole which or- 
ders at low temperatures in Am02- 

Experimentally it has been considered that the CEF 
ground state of AmC>2 is iy ground state. However, we 
believe that it is still meaningful to examine the exper- 
imental results on the basis of the T% quartet ground 
state, although it may be difficult to perform the micro- 
scopic experiments of Am0 2 . 

Finally, let us provide a comment on the simplifica- 
tion of the model. In this paper, since we have consid- 
ered only single conduction band, there exists residual 
entropy in the results. In actuality, it should be finally 
released when we consider t\ u and t2 U conduction bands. 
This point is also related to the relevant multipole mo- 
ment when we consider the ordered state in the periodic 
systems. In this sense, the present results are qualita- 
tive, but they include the actual multipole which forms 
ordered state. 

In summary, we have discussed the multipole state of 



actinide dioxides due to the evaluation of the multipole 
susceptibility of the Anderson model. When Coulomb 
interaction and spin-orbit coupling have been appropri- 
ately changed, it has been found that multipoles includ- 
ing quadrupole and octupole are relevant to Am02- It 
is believed that multipole ordering can be detected in 
Am02 in future experiments. 
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